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Electron localization is the tendency of an electron in a many-body system to exclude other 
electrons from its vicinity. Using a new natural measure of localization based on the exact many- 
electron wavefunction, we find that localization can vary considerably between different gronnd-state 
systems, and can also be strongly disrupted, as a function of time, when a system is driven by an 
applied electric field. We use our new measure to assess the well-known electron localization function 
(ELF), both in its approximate single-particle form (often applied within density-functional theory) 
and its full many-particle form. The full ELF always gives an excellent description of localization, but 
the approximate ELF fails in time-dependent situations, even when the exact Kohn-Sham orbitals 
are employed. 

PACS numbers: 72.15.Rn, 71.15.Mb, 31.15.A- 


Density functional theory (DFT) [1] replaces the many- 
body (MB) wavefunction with the electron density as 
its fundamental variable. This transformation, however, 
provides new challenges as, in the non-interacting world 
of Kohn-Sham (KS) electrons, the many-body behavior 
of electrons is concealed within the exchange-correlation 
(xc) potential. Electron localization, describing the po¬ 
sition entanglement of electrons, is one such property. 
Pauli exclusion tends to localize like-spin electrons in sep¬ 
arate regions. The Coulomb repulsion further increases 
the tendency to localize. Localization, however, increases 
the kinetic energy, and the balance between these fac¬ 
tors makes it challenging to quantify electron localiza¬ 
tion. An understanding of electron localization is useful 
chemically, linking with the concepts of chemical bonds 
and localized electron pairs [2], and also describes a fun¬ 
damental aspect of electron correlation that approximate 
DFT functionals should take into account [3H5]. 

We stress that localization is a true many-body prop¬ 
erty of electrons, dependent on the positions of all elec¬ 
trons, and not accessible through the spatial character of 
any one KS orbital. 

In this Rapid Communication, we explore the varia¬ 
tion in localization for a range of ground-state and time- 
dependent systems. Using a new measure of localization, 
calculated from the MB wavefunction, we find a supris¬ 
ing degree of delocalization when we perturb the system 
with a strong electric held. We then compare it with 
the usual formulation of the electron localization function 
(ELF), often used to characterize electron localization in 
DFT calculations [6]. We hnd that this ELF correctly 
describes localization across a range of ground-state sys¬ 
tems, but fails in time-dependent problems where delo¬ 
calization occurs due to electron collisions. These errors 
are not present, however, if the ELF is formulated from 
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the true pair density. 

To calculate our measure we partition the system into 
distinct regions, each containing charge corresponding to 
one electron, and examine the actual distribution of the 
N electrons among these regions [7j. The probability 
that each region contains only one electron is readily cal¬ 
culated as 

p = J M{ri,...,rN) |T'(ri,... ,rAr)|^ d^ri... d^rAT, 

( 1 ) 

where the mask function, M, dehned as 1 if all electrons 
are in different spatial regions and otherwise 0, interro¬ 
gates the occupation of these regions. 

To form our regional electron localization measure 
(RELM), we scale p with reference to the probability 
Po = Nl/N^ of finding exactly one electron in each re¬ 
gion for an ideal delocalized and uncorrelated state: 

RELM=^^^^. (2) 

1 - To 

This makes our measure 0 when a system is ideally delo¬ 
calized, and 1 when fully localized. As the probability of 
simultaneously finding only one electron in each region 
increases, so does our measure of localization. 

Unfortunately, in most cases the MB wavefunction is 
too expensive to calculate. The traditional method of 
approximating the localization in a system is the ELE 
of Becke and Edgecombe [8]. Originally developed for 
Hartree-Fock (HF) calculations, the method can also be 
applied to Kohn-Sham orbitals. The ELF is based on the 
quantity Da-, defined by the Taylor series 

pK.nAr,s)=^^^^+Ois^), (3) 

where Pcond ^ conditional probability, the probability 
given that an electron has been found at position r that 
a second electron will be found at distance s from this 
position, d is the spatial dimensionality of the system 
and (T is a spin index. Da characterizes the probability 
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of finding a second electron very close to the reference 
electron and is a measure of localization in its own right. 

Later work by Dobson [5] provides an equivalent def¬ 
inition of D„ (r) to that given in Eq. allowing this 
quantity to be calculated directly from the pair density, 
using the equation 


Da (r) 


V?,n2(r,r') 

2n(r) 


(4) 


where n 2 {r,r') is the pair density. Using this expres¬ 
sion, it is possible to calculate the ELF directly from the 
wavefunction; we term this the “exact ELF”. 

As has a strong dependence on the local density, 
it is not directly interpreted easily. To produce the ELF, 
Becke and Edgecombe scaled Da as 



This expression compares the local value of Da with that 
of a homogeneous electron gas (HEG) of the same lo¬ 
cal density. Da,a- Hence, ELF ranges from zero to one, 
where one represents total localization and a half rep¬ 
resents the degree of localization in a HEG of the same 
density. 

In general, the pair density is unknown, so the exact 
definition of Eq. can not be used. Becke and Edge¬ 
combe found an approximate expression to calculate Da 
in terms of single-particle orbitals (p: 


Da 




( 6 ) 


where na is the electronic density. We use Hartree atomic 
units (a.u.) here and henceforth. This equation is exact 
in a HE treatment, and is also commonly applied to KS 
orbitals [101 • We term this the “approximate ELF”. 

As shown, Becke and Edgecombe’s approach relies on 
two main assumptions: that Da (r) is an effective local 
description of localization, and that the approximation 
to it given by Eq. is satisfactory. If either of these as¬ 
sumptions fail, the ELF will give a misleading picture of 
electronic behavior. ELF calculations have been widely 
used, but little is known about their accuracy. 

ELF is a local measure of localization and in order 
to look at the localization of systems as a whole we de¬ 
fine the average ELF, weighted according to the density, 
(ELF) = ^ J ELF (r) n (r) , which can range from 

zero to one. 

All the results presented in this Letter were calculated 
using the iDEA code m- The Schrodinger equation is 
solved exactly for ID systems of spinless electrons to find 
the MB wavefunction, both for ground-state and for time 
propagation after an electric field is applied. This pro¬ 
vides the charge density for all our model systems. Next, 


the exact KS orbitals that reproduce the density are cal¬ 
culated m- As is standard in ID, we use a softened 
Goulomb interaction. 

Some components of ELF are dimensionally depen¬ 
dent. The appropriate changes in 2D have already been 
discussed [T^ and we follow this approach. is di¬ 
mensionally dependent and we use Da H = in 

this work. Previously ID ELF has been applied to cou¬ 
pled electron-nuclear wavefucntions |13j . The partitions 
needed for ID HELM are naturally chosen, although it is 
important to update them in a time-dependent system. 

First we look at a family of three-electron double-well 
potentials and calculate the localization of their ground- 
states. These external potentials are defined as 

V (x) = ax^^ — /3x‘^. (7) 

A constant value of a = 5 x 10“^^ a.u. was used, while 
the value of /3 was varied: when /3 = 0 there is no barrier 
in the potential, and as it is increased the height of the 
barrier grows. Three examples of the potentials and den¬ 
sities for different values of /3 are shown in Fig. This 
family of wells is interesting as the double-well potential 
only provides two natural sites for the three electrons to 
occupy. 

Fig- 0^) shows some of the ground-state potentials 
that make up this family and the effects of the poten¬ 
tial barrier on the electrons. At first, the middle electron 
stays in the barrier region owing to the Coulomb repul¬ 
sion, as its height increases in (i) to (ii). As the pres¬ 
ence of the barrier disperses the central electron across 
the system, it also drives the outer electrons toward the 
boundaries of the system, acting to increase the localiza¬ 
tion of the system rather than decrease it (see Fig. 0b)). 
As the strength of the barrier is increased further, it be¬ 
comes energetically favorable for the central electron to 
move into the two side wells, as in (iii), reducing local¬ 
ization. 

In Fig. 0b), the three measures agree on how the local¬ 
ization of this family of potentials varies. The similarity 
between RELM and the exact average ELF is striking 
as the two methods are based on different mathematical 
interpretations of localization: RELM is scaled proba¬ 
bilistically and ELF is scaled with respect to the HEG 
as reference. Approximating Da does lead to a system¬ 
atic lowering of the calculated (ELF), but still yields the 
correct trend across the range of localizations calculated. 

We next investigate the information contained in Da 
in more detail. Fig. shows an example plot for a two 
electron system (specified in the caption), demonstrating 
the approximation made in Eq. in practice. The con¬ 
ditional probability Pcond is plotted for a chosen value 
of r. As demanded by the Pauli exclusion principle, the 
probability of finding a second electron at the same po¬ 
sition (s = 0) is zero. The function then shows a peak at 
the most likely separation the second electron is found at 
(we would normally expect as many peaks as there are 
remaining electrons in the system). 
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FIG. 2. (color online) A plot of electron separation s against 
the conditional probability PcL^d {xi,s) where xi is fixed as 
the density maximum (2.64 a.u.) for the two electron ground- 
state of the potential of Fig. [^a)(i). The conditional proba¬ 
bility (blue solid) is plotted with the ELF level approximation 
to it (green dashed). Insert: magnification of short s behav¬ 
ior. Da does not contain any long range information and only 
correctly characterizes short distances, in this example only 
~ 0.6 a.u. Our strong ELF results suggest that this neglected 
long s behavior is not important. 


FIG. 1. (Golor online) 3-electron double wells — (a) Plots of 
the external potentials (dashed blue) of three selected wells as 
the barrier height is increased. The ground-state charge den¬ 
sities (solid green) of these potentials are shown. The local¬ 
ization regions used in the RELM calculations are also shown 
(dotted black); each contains exactly one electron’s worth of 
charge, (b) The localization of the family of potentials, cal¬ 
culated using the three methods introduced - RELM, exact 
average ELF and approximate average ELF. All three mea¬ 
sures agree how the barrier influences localization. Triangles 
indicate the values of /3 for which the potentials are plotted 
in (a). 


Also shown on the plot is the approximated version of 
this conditional probability that is used in ELF calcu¬ 
lations. As shown in Eq. the ELF approximates this 
conditional probability as (r) /d and this is shown 
in the plot where has been calculated using Eq. As 
shown in the inset in Fig. this approximation is only 
effective over very short electron separations. RELM cal¬ 
culations use information over all s, so the agreement be¬ 
tween our ELF and RELM calculations suggests that this 
longer-range s behavior is not an important ingredient for 
a localization measure to contain. 

Next we look at a time-dependent system. In the time- 
dependent regime the approximate ELF is modified by 
the addition of an extra term to Eq. producing the 
time-dependent ELF (TDELF) [T3]. This equation be¬ 
comes 


Da 


Na 


|VCI'- 


1 iVn^l 


3(7 

Tt(j 


( 8 ) 


where ja is the current density. 

For time-dependence, we return to the potential well 
having /3 = 0, shown in Fig. [^a)(i), and this time place 
two electrons in it. As before, we find the ground-state 
wavefunction of the well. Then at time < = 0 we apply 
a perturbing uniform d.c. electric field, driving the elec¬ 
trons strongly towards the right-hand well, causing them 
to “collide”. We also look at the non-interacting system 
with Pauli exclusion but no Coulomb interaction. 

Fig. i (a) shows how the localization of the interact¬ 
ing electrons changes along the 80 a.u. of the simulation. 
Broadly, these results show significant changes in local¬ 
ization over time. This is in contrast to the notion that 
localization is a persistent characteristic of a system. 

On both plots the expectation value of the electron 
separation (s) is shown. The large localization drops oc¬ 
cur around the two electrons becoming closer together. 
These rapid drops in localization also appear in other 
systems we have investigated and seem to be a common 
feature of electron collisions. For our systems, (s) seems 
to be a fair indicator of localization in its own right. 

Fig. §b) shows the calculation repeated with no 
Coulomb interaction. The behaviour is broadly similar, 
but interaction seems to exaggerate some features and 
suppresses others. This comparison strengthens the ar¬ 
gument that Pauli exclusion is the main driver of local¬ 
ization. 

Again, for both plots RELM and exact ELF show the 
same trends, although both measures show some unique 
features. The close agreement between the two measures 
still holds in a time-dependent context. We also note 
the agreement between exact ELF and RELM is reduced 
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FIG. 3. (color online) Localization measures with time for 
the perturbed well, (a) The interacting system. RELM 
(solid blue) and exact average ELF (dashed green) continue to 
show good agreement. The expectation value of the electron 
separation (long dashed red) is also shown on the second axis 
and shows similar features, (b) The same plot for the non¬ 
interacting system. The behaviour of the system is similar, 
showing that Pauli exclusion is the main driver of localization. 
RELM and ELF differ more signihcantly 


FIG. 4. (color online) The curves in Fig. [^for a shorter time 
interval, with the addition of approximate ELF calculations. 
For the interacting system in (a), the average approximate 
TDELF (dotted-dashed purple) and the GS approximate ELF 
(dotted orange) both delocalize the electrons too soon. The 
TDELF, though an improvement on the original approxima¬ 
tion, has a large spurious drop in localization around t = 5 
a.u. As shown in (b), this approximation performs better 
without interaction as the problematic drop is weaker. 


when the Coulomb interaction is turned off. 

Fig. a shows approximate ELF calculations for the 
same systems: both the original formulation and the 
time-dependent extension. (These are shown for a 
shorter time interval due to the numerical challenge in 
calculating exact KS orbitals for later timesteps.) For 
the interacting system, the original (time-independent) 
ELF formulation performs very poorly and for most of 
the calculation shows the system as erroneously delocal¬ 
ized. The extra current term in the TDELF makes a 
significant improvement. At around 5 a.u., however, it 
too shows an unphysical drop in localization which anal¬ 
ysis shows to occur in the left-hand well, predominantly 
associated with a large increase of the first term of Eq. 
in the left half of the system which is not compensated for 
by the other terms. Instead, the real system delocalizes 
both electrons at 12 a.u., later in the simulation. Clearly, 
the TDELF’s ignorance of correlation (beyond Pauli ex¬ 
clusion) in the wavefunction is limiting its description of 
localization. 


The approximate TDELF performs better when there 
is no electron-electron interaction [Fig. I^b)]. It still 
shows a slow drop when the localization is staying con¬ 
stant, but this erroneous drop is significantly weaker. 

One failure of the present approximation, concealed in 
the definition of ELF, is that it is not positive definite, 
leading to non-physical negative values of Do-. If these 
values are set to zero, a small improvement in accuracy 
is achieved. Negative values should serve as a warning 
that the method is not performing reliably m- 

In conclusion, we have studied electron localization 
using a variety of measures accross a range of ground- 
state and time dependent systems. Our results show 
the strength of the ELF approach, despite its focus on 
short range exclusion. We further find that the usual 
approximate ELF provides good results for our ground- 
state systems, notwithstanding its simplicity and neglect 
of correlation. In contrast, time-dependent systems can 
often become surprisingly delocalized as electrons collide 
with one another, and in this case the simple approxi- 
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mate ELF is no longer adequate. 
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